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We investigate the propagation of electron wavepackets in single-walled carbon nanotubes via a Lindblad- 
based density-matrix approach that enables us to account for both dissipation and decoherence effects induced 
by various phonon modes. We show that, while in semiconducting nanotubes the wavepacket experiences the 
typical dispersion of conventional materials, in metallic nanotubes its shape remains essentially unaltered, 
even in the presence of the electron-phonon coupling, up to micron distances at room temperature. 
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One of the most fascinating challenges in today’s nan¬ 
otechnology is to control the dynamics of electron waves 
with an accuracy comparable to the one available in 
optics. 1 To this purpose, it is crucial to generate se¬ 
quences of wavepackets that propagate coherently with¬ 
out overlapping to each other. Although single-electron 
pumps have been implemented, 2-4 electron waves exhibit 
a major limitation as compared to their electromagnetic 
counterpart: While the photon velocity is independent 
of its wavevector, the electron group velocity in conven¬ 
tional materials (with a parabolic-band structure) de¬ 
pends on it, causing the well known intrinsic disper¬ 
sion of the electron wavepacket, even in the absence of 
scattering processes. Thus, despite various proposals, 5-8 
the realization of single-electron wavepackets of control¬ 
lable shape and phase that propagate ballistically in low¬ 
dimensional conductors still remains an open problem. 

In metallic single-walled carbon nanotubes (SWNTs), 
in graphene and in the recently discovered topological in¬ 
sulators, however, electrons behave as massless relativis¬ 
tic fermions and, just like photons, exhibit a linear spec¬ 
trum, with the Fermi velocity vf ~ 10 6 m/s replacing 
the speed of light c. This property makes such materials 
ideal candidates for an electronic alternative to photon- 
based quantum information processing. In graphene, for 
instance, electron supercollimation has been predicted to 
occur when an external static and long-range disorder 
is suitably applied. 9,10 SWNTs, nowadays synthesized 
with high accuracy 11-13 and used as versatile building 
blocks of high-performance electronic devices, 14-16 are 
even more promising. 

Although in principle an electron wavepacket can prop¬ 
agate along a metallic SWNT maintaining its initial 
shape, any implementation in realistic devices requires to 
determine whether and to what extent such property is 
affected by scattering processes. The high level of purity 
reached in carbon-nanotube synthesis 12 enables one to 
fairly neglect scattering with vacancies or defects. Fur¬ 
thermore, charged-impurity potentials arising from the 
substrate are typically long-ranged over the C-C distance, 
and are well known to produce no backscattering, due to 


the Klein tunneling effect related to the two-fold sub¬ 
lattice degree of freedom characterizing C-based mate¬ 
rials. At low temperatures, electron-electron scattering 
has been shown to lead to Coulomb-blockade 17 and to 
Luttinger-liquid behaviors. 18 However, at intermediate 
and room temperature, experiments indicate that the 
main source of scattering arises from electron-phonon 
coupling, 19-21 whose effects have attracted the interest 
of theorists in last years. On the one hand, models based 
on a classical-like treatment of the electron-phonon cou¬ 
pling (as an external oscillating potential) 22,23 enable one 
to obtain the time-dependent evolution of single wave 
packets and to extract the linear conductance by a suit¬ 
able average over the initial state. These approaches, 
however, fail in capturing the intrinsically dissipative na¬ 
ture of the phonon bath. On the other hand, treat¬ 
ing electron-phonon coupling in SWNTs via Boltzmann- 
equation schemes 24 does not allow one to account for 
electronic phase coherence. 

In this Letter we investigate the impact of electron- 
phonon coupling on the dynamics of an electron 
wavepacket propagating in SWNTs; to this end, we em¬ 
ploy a recently developed density-matrix approach that 
enables us to account for both energy dissipation and 
decoherence effects. We demonstrate that, while in semi¬ 
conducting SWNTs an electronic wavepacket spreads 
already for a scattering-free propagation, in metallic 
SWNTs the shape of the wavepacket is essentially un¬ 
altered, even accounting for electron-phonon couplings. 
Differently from predictions based on the conventional 
relaxation-time approximation (RTA), our results thus 
indicate that metallic SWNTs can realistically be utilized 
as an electron-based platform for information transfer. 

The low-energy carrier Hamiltonian of a SWNT de¬ 
couples into two valleys around the K and K 7 points, de¬ 
scribed by Dirac matrices 77k = Hvfct ■ k and H k' = 
— Uvfct* • k, respectively, where er = (, a x ,a y ) denote 
Pauli matrices acting on the twofold sublattice space, 
and k = (k^ u ,k) the carrier wavevector. 25 Here k de¬ 
notes the continuous component along the SWNT axis 
(||), whereas h^ v = (n + w/3)/R is the discrete com- 
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ponent along the circumference (_L), with n denoting the 
electron subband, v = ±1 for the K/K' valley, R the nan¬ 
otube radius, and the index v = 0, ±1 is defined through 
the relation exp[zK • C] = exp[—*K' • C] = exp[27rw/3], 
where C is the vector rolling the graphene lattice into the 
SWNT. Since typical subband energy separations are of 
the order of eV, we shall focus on the lowest energy sub¬ 
band (n = 0), whose energy spectrum is given by 


e a = bhvF\/k 2 + (is/3R) 2 , (1) 


where a = (fc, b, v, u) is the quantum-number multi-label 
and b = ±1 for the conduction and valence band, respec¬ 
tively. The related electron eigenfunctions are <j> a { r) = 
(4tt RL)- 1 / 2 (e- iv6k / 2 ,bve +ivekl2 )e ik - r , where L denotes 
the nanotube length and 6 k the polar angle of k. While 
for v ^ 0 the energy spectrum is gapped (semiconducting 
nanotube) and near k = 0 is parabolic-like similarly to 
conventional semiconductors, for v = 0 the spectrum is 
gapless (metallic case), and the typical massless Dirac- 
cone structure is recovered. All armchair and (3n,0) 
zigzag SWNTs are remarkable examples of the metallic 
case. 25 

The SWNT phonon spectrum near the K and K' 
points includes 25,26 i) longitudinal(L) stretching modes, 
characterised by an acoustic(A) branch u.> q ,LA = vs|g| 
with vs — 1.9 • 10 4 m/s and an optical(O) branch with 
huiLO — 0.2 eV; ii) breathing(Br) modes orthogonal 
to the nanotube surface, with a roughly g-independent 
spectrum hujBr — 0.14eVA/i?; (iii) transverse(T) twist¬ 
ing modes, characterised by an acoustic branch with 
fjJq,TA = vt\q\ with vt — 1.5 • 10 4 m/s and an optical(O) 
branch with Tiloto — 0-2 eV. For each mode the electron- 
phonon coupling is described by a 2 x 2 matrix acting on 
the sublattice space, which exhibits both diagonal (intra¬ 
sublattice) and off-diagonal (inter-sublattice) entries for 
the LA, TA and Br modes, 26 and is purely off-diagonal 
for optical modes. 27 We neglect zone-edge phonon modes 
that lead to inter-valley scattering because they may be¬ 
come important only at very high energies. 

In order to account for energy dissipation and deco¬ 
herence induced by the nanotube phonon bath on the 
otherwise phase-preserving electron dynamics, we have 
applied to the carbon nanotube the formalism introduced 
in Ref. [28] via a numerical solution of the Lindblad-based 
non-linear density-matrix equation (LBE) 
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Here, the first term describes the scattering-free prop¬ 
agation, whereas the second term is a non-linear scat¬ 
tering superoperator expressed via generalized scatter¬ 
ing rates V^ ia2 a i a i , whose explicit form is microscopi¬ 


cally derived from the electron-phonon Hamiltonians de¬ 
scribed above, with s labelling the various phonon modes. 
The fully quantum-mechanical density-matrix equation 
(2) enables us to go beyond the conventional Boltzmann 
transport equation, which is recovered in the diagonal 
limit (p a ia 2 = f ai S aia2 ), 29 where the generalized scat¬ 
tering rates reduce to the semiclassical rates provided by 
the standard Fermi’s golden rule, P b a , = V b Qa a , a ,. 

In order to show that carbon nanotubes can be utilised 
as quantum-mechanical channels for the non-dispersive 
transmission of electronic wavepackets, we have per¬ 
formed simulated experiments where the shape of an ini¬ 
tially prepared wave packet is monitored while it evolves 
under the effect of the phonon bath. Any initial state 
can always be written as p = p° + 5p, where p° is the 
homogeneous equilibrium state and Sp describes a lo¬ 
calised excitation, whose spatial shape (e.g. Gaussian) 
can in principle be generated experimentally via a prop¬ 
erly tailored optical excitation. While the description of 
the specific optical generation is beyond the aim of the 
present paper, the localisation of the initial wave packet 
is the crucial aspect in our analysis. A simple choice 
that captures this essential physical ingredient is an ini¬ 
tial state described by an intra-valley conduction-band 
density matrix p aia2 = S Vl , V2 S bl) iS b2t ip kl k 2 , where 


= • ( 3 ) 

Here /£ is the Fermi-Dirac distribution of the conduction- 
band states fc, and the parameter £ plays the role of a 
delocalization length. Indeed for £ —> oo the spatially 
homogeneous equilibrium state p% k = f kl dk 1 k 2 is 
recovered, whereas for finite £ the excitation Sp consists 
of interstate phase coherence (intraband polarization) 
determined by the parameter l. In particular, in 
the limit £ —> 0 one obtains the maximally localized 
wavepacket. Moreover, a distinguished feature of the 
initial condition in (3) is the absence of nonequilibrium 
diagonal contributions ( p kk = /£), which implies that 
energy dissipation and decoherence will affect the non¬ 
diagonal contributions only; this is the typical situation 
produced by a weak interband optical excitation. 

We start our discussion analysing the scattering-free 
propagation of the initial state (3), e.g. switching off the 
electron-phonon interaction in Eq. (2). Figure 1 shows 
the spatial carrier density n(r ||) at different times for the 
initial state (3), taken at room temperature and in the 
maximally-localized limit (£ —> 0). The cases of semicon¬ 
ducting and metallic nanotube are shown in the upper 
and lower panels, respectively. For the semiconducting 
nanotube {v ^ 0) the dispersion relation (see Eq. (1)) at 
small k reduces to the parabolic spectrum of conventional 
semiconductor materials and gives rise to the typical 
classical-like diffusion process, preventing any informa¬ 
tion transfer via electronic wavepackets. In contrast, for 
the case of the metallic nanotube (v = 0), characterised 
by a linear dispersion, the initial charge peak splits into 
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FIG. 1. (Color online) Room-temperature scattering-free evo¬ 
lution in a semiconducting (10,0) SWNT (u = 1) (upper 
panel) and a metallic (12,0) (u = 0) SWNT (lower panel): 
Charge density of the maximally-localized (£ —¥ 0) electronic 
wavepacket in Eq. (3) as a function of the position along the 
SWNT axis at three different times: t = 0 fs (solid), t = 50 fs 
(dashed) and t = 100 fs (dashed-dotted). 


two identical and shape-preserving components which 
travel in opposite directions with velocity ±vf (lower 
panel), i.e., n{r y, t) = n + (r || — vpi) + n~{r || + upf). This 
is the non-dispersive propagation scenario we are looking 
for. 

We thus focus on the metallic nanotube and switch on 
the electron-phonon coupling to address the crucial ques¬ 
tion of whether and how energy dissipation and decoher¬ 
ence modify such ideal dispersion-free scenario. To this 
end, we have performed a set of simulated experiments 
based on the LBE (2), including LA, TA, Br, LO and TO 
phonons. In order to expand the space scale with respect 
to the ideal scenario displayed in Fig. 1, here the delocal¬ 
ization length i in the initial condition (3) is chosen such 
to get a FWHM value of the initial peak of 0.4/rm. 

Figure 2(a) shows the wavepacket propagation 
in a (12,0) SWNT at room-temperature. The electron- 
phonon coupling, accounted for by our LBE simulation 
(solid curve), does not significantly alter the shape with 
respect to the ideal scattering-free result (dashed curve), 
so that the transmission is essentially dispersionless up 
to the micrometric scale. In order to understand the 
origin of such shape-preserving dynamics, we first ob¬ 
serve that in the electron-phonon coupling a natural 
distinction arises between forward and backward scat¬ 
tering processes, where the initial and final electronic 
states have the same and opposite velocity sign, respec¬ 
tively. A semiclassical argument based on Boltzmann 
theory leads to expect that in metallic SWNTs, where 



position (pm) 



FIG. 2. (Color online) (a) Room-temperature charge-density 
evolution of an electronic wavepacket in a metallic (12,0) 
SWNT as a function of the position r\\ along the SWNT axis, 
at three different times: t = 0 ps (left peaks), t = 1 ps (central 
peaks) and t = 2 ps (right peaks). The solid lines show the 
effects of the carrier-phonon coupling, accounted for by the 
LBE (2), compared to the scattering-free propagation (dashed 
lines). The RTA results (dot-dashed lines) are also plotted for 
comparison. Here the wavepacket delocalization length i in 
Eq. (3) is chosen such that the FWHM initial value is 0.4 
pm. (b) Room-temperature total scattering rates in a metal¬ 
lic (12,0) SWNT as a function of the axial energy = hvF\k\ 
for various electron-phonon mechanisms s =LA (dotted line), 
TA (short-dotted line), Br (dot-dot-dashed line), LO (short- 
dot-dashed line), TO (dot-dashed line), as well as for their 
sum (solid line) and its backward component (dashed line). 
Note that the s =LO and s =TO lines almost coincide. 


the group velocity can only take two fc-independent val¬ 
ues ±wf, forward processes never contribute to the dis¬ 
persion of wavepackets, which is thus due to backward 
processes only. In our fully quantum-mechanical LBE 
approach the situation is less trivial, as the semiclassi¬ 
cal transition k —> k + q between purely diagonal en¬ 
tries of the density matrix extends to off-diagonal quan¬ 
tum transitions (^ 1 ,^ 2 ) —> (fci + q, &2 + <?), so that for¬ 
ward and backward processes can in principle interplay 
in Eq. (2). In semiconducting materials such transitions 
can lead to scattering nonlocality and quantum diffusion 
speed-up phenomena; 29 moreover, in Luttinger liquids 
the forward component of the electron-phonon coupling 
can lead to Wentsel-Bardeen instabilities of the electron 
propagator. 30 However, for a metallic SWNT, a detailed 
investigation (not reported here) shows that, due to a 
strong cancellation between in- and out-scattering terms, 
the forward-process contribution to the scattering super¬ 
operator in (2) has an extremely small impact on the 
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wavepacket propagation. This enables us to support on 
firmer grounds the conclusion that in metallic SWNTs 
the wavepacket dispersion (solid line in Fig. 2(a)) orig¬ 
inates mainly from backward processes. At room tem¬ 
perature, the latter are ascribed to breathing phonons 
and to different acoustic modes, depending on the type 
of SWNT: for a zigzag SWNT, like the one considered 
in Fig.2(a), they are due to LA modes, whereas for an 
armchair SWNT they are due to TA modes. 

Notably, a strong suppression of the wave packet would 
be predicted by a conventional RTA model, as shown by 
the dash-dotted line in Fig. 2(a). Indeed, such RTA is 
formulated in terms of total scattering rates T* (from 
an initial state a to any final state) for each phonon 
mode s. The energy dependence of T* , displayed by the 
solid line in Fig.2(b), shows that the most relevant scat¬ 
tering mechanisms are due to LA and Br modes, lead¬ 
ing to the wrong conclusion that the typical phonon- 
induced wavepacket dispersion time-scale is about 2 ps, 
in accordance to the RTA result (dash-dotted line in 
Fig. 2(a)). However, the dominant contributions of LA 
and Br modes originate from forward processes. Indeed, 
in the 2x2 sublattice matrix characterising the SWNT 
electron-phonon interaction for both acoustic and breath¬ 
ing modes, the coupling constant of the diagonal en¬ 
tries is known to dominate by a factor ~ 15 — 25 the 
off-diagonal ones. 25,26 Because backward scattering pro¬ 
cesses can only arise from off-diagonal entries, the total 
backward component of the scattering rate is at least 
one order of magnitude smaller than the total rate (see 
dashed line in Fig. 2(b)). For the reasons discussed 
above, namely a strong cancellation between in- and out- 
scattering terms (absent within the RTA model), these 
processes have in fact negligible impact on the diffusion 
process. We stress that also LO modes induce backward 
processes; however, due to their high phonon energy, 19 
their effect is negligible. We also notice that, at lower 
temperatures, the impact of electron-phonon coupling de¬ 
creases and the diffusion time becomes even longer (for 
instance, at 77 K the wavepacket travels for more than 
one micron with a signal attenuation less than 1%). 

Finally, we emphasise that, compared to wavepacket 
propagation in parabolic-like systems (e.g., conventional 
semiconductor materials 29 or semiconducting nanotubes 
(see upper panel in Fig. 1)), in a metallic SWNT the im¬ 
pact of phonon-induced diffusion on the wavepacket prop¬ 
agation is very limited; indeed, while for a parabolic-like 
system the free-diffusion spatial-broadening displayed in 
the upper panel of Fig. 1 is further increased due to its 
coupling with various dissipation channels, here the lin¬ 
ear nature of the spectrum mostly preserves the shape 
of the initial peak, which is thus free to move at Fermi- 
velocity remaining almost unaltered despite the presence 
of dissipation/decoherence processes. Our results thus 
indicate that metallic SWNTs are a promising platform 
to realise quantum channels for the non-dispersive trans¬ 


mission of electronic wavepackets. 
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